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We review the known properties of molecular outflows from low- and high-mass 
young stars. General trends among outflows are identified, and the most recent studies 
on the morphology, kinematics, energetics, and evolution of molecular outflows are 
discussed, focusing on results from high-resolution millimeter observations. We 
review the existing four broad classes of outflow models and compare numerical 
simulations with the observational data. A single class of models cannot explain flie 
range of morphological and kinematic properties that are observed, and we propose 
a possible solution. The impact of outflows on their cloud is examined, and we 
review how outflows can disrupt their surrounding environment, through the clearing 
of gas and the injection of momentum and energy onto the gas at distances from 
their powering sources from about 0.01 to a few pc. We also discuss the effects of 
shock-induced chemical processes on the ambient medium, and how these processes 
may act as a chemical clock to date outflows. Lastly, future outflow research with 
existing and planned millimeter and submillimeter instruments is presented. 



1. INTRODUCTION 

As a star forms by gravitational infall, it energetically 
expels mass in a bipolar jet. There is strong evidence for a 
physical link between inflow and outflow and that magnetic 
stresses in the circumstellar disk-protostar system initially 
launch the outflowing material (see chapters by Pudritz et 
ai; Ray et al. \ and Shang et al). The ejected matter can ac- 
celerate entrained gas to velocities greater than those of the 
cloud, thereby creating a molecular outflow. Outflows can 
induce changes in the chemical composition of their host 
cloud and may even contribute to the decline of the infall 
process by clearing out dense gas suiTounding the protostar. 



In addition, molecular outflows can be useful tools for un- 
derstanding the underlying formation process of stars of all 
masses, as they provide a record of the mass-loss history of 
the system. 

Protostellar outflows can be observed over a broad range 
of wavelengths, from the ultraviolet to the radio. In this 
review we will concentrate on the general characteristics 
and properties of molecular outflows, the entrainment pro- 
cess, and the chemical and physical impact of outflows on 
the cloud that are mainly detected through observations of 
molecular rotational line transitions at millimeter and sub- 
millimeter wavelengths. At these wavelengths the obser- 
vations mainly trace the cloud gas that has been swept-up 
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by the underlying protostellar wind, and provide a time- 
integrated view of the protostar's mass-loss process and its 
interaction with the surrounding medium. 

2. GENERAL OUTFLOW PROPERTIES 

Over the last 10 years, millimeter interferometers have 
allowed the observation of molecular outflows at high angu- 
lar resolutions (~ 1 to 4"), while the capability to observe 
mosaics of several adjacent fields has enabled mapping of 
complete outflows at those resolutions. Such interferomet- 
ric observations give access to the internal structure of the 
gas surrounding protostars, and can disentangle the mor- 
phology and dynamics of the different elements that are 
present (i.e., protostellar condensation, infalUng and out- 
flowing gas). These high resolution observations have been 
critical to the discovery of the kinematics and morphology 
of outflows from massive OB (proto)stars, which are typi- 
cally more than a kiloparsec away. 

General trends have been identified in molecular out- 
flows from both low- and high-mass protostars, even though 
they display a broad diversity of sizes and shapes. These 
properties have been identified mostly using single-dish 
and interferometer observations of the CO Unes. Molec- 
ular outflows exhibit a mass-velocity relation with a bro- 
ken power law appearance, dM{v)/dv oc v~'', with the 
slope, 7, typically ranging from 1 to 3 at low outflow ve- 
locities, and a steeper slope at higher velocities — with 7 
as large as 10 in some cases (e.g., Rodriguez et al, 1982; 
Lada andFich, 1996; Ridge and Moore, 2001). The slope 
of the mass-velocity relation steepens with age and energy 
in the flow {Richer et al, 2000). The velocity at which 
the slope changes is typically between 6 and 12 kms~^ 
although outflows can have CO break velocities as low as 
about 2 km s~^ and, in the youngest CO outflows, it can be 
high as 30 km s~^ (see, e.g.. Richer et al., 2000, and refer- 
ences therein). The mass, force, and mechanical luminosity 
of molecular outflows correlate with bolometric luminos- 
ity {Bally and Lada, 1983; Cabrit and Bertout, 1992; Wu 
et al, 2004), and many fairly collimated outflows show a 
linear velocity-distance relation, typically referred to as the 
"Hubble-law", where the maximum radial velocity is pro- 
portional to position (e.g., Lada and Fich, 1996). Also, the 
degree of coUimation of outflows from low- and high-mass 
systems appears to decrease as the powering source evolves 
(see below). 

These observed general trends are consistent with a com- 
mon outflow/infall mechanism for forming stars with a wide 
range of masses, from low-mass protostars up to early B 
protostars. Although there is evidence that the energetics 
for at least some early-B stars may differ from their low- 
mass counterparts, the dynamics are still governed by the 
presence of hnked accretion and outflow. A few young O 
stars show evidence for accretion as well although this is 
not as well established as for early-B stars (e.g., van der 
Tak and Menten, 2005; chapter by Cesaroni et al.). 



2.1 Outflows from low-mass protostars 

Since their discovery in the early eighties, molecular out- 
flows driven by young low-mass protostars (i.e. typically 

< 1 M0) have been extensively studied, giving rise to a 
detailed picture of these objects (see, e.g., the reviews by 
Richer et ah, 2000; Bachiller and Tafalla, 1999, and refer- 
ences therein). The flows typically extend over 0.1-1 par- 
sec, with outflowing velocities of 10-100 kms~^. Typi- 
cal momentum rates of 10~^ M© km s~^yr~^ are observed, 
while the molecular outflow mass flux can be as high as 
10~^ M0yr~^ {Bontemps et al., 1996). Particular interest 
has been devoted to the outflows driven by the youngest, 
embedded protostars (age of a few 10'^ to a few 10* years, 
the Class objects). These sources are still in their main 
accretion phase and are therefore at the origin of very pow- 
erful ejections of matter 

2.7.7. Molecular jets. The colhmation factor (i.e., 
length/width, or major/minor radius) of the CO outflows, 
as derived from single-dish studies, range from ^3 to 
>20. There is however a clear trend of higher colhma- 
tion at higher outflowing velocities (see, e.g., Bachiller 
and Tafalla, 1999). Interferometric maps have revealed 
even higher colhmation factors, and, in some cases, high- 
velocity structures that are so coUimated (opening angles 

< a few degrees) that they can be described as "molecular 
jets". 

HH211 is the best example to date of such a molecu- 
lar jet {Gueth and Guilloteau, 1999). At high-velocity, the 
CO emission is tracing a highly-coUimated hnear structure 
that is emanating from the central protostar. This CO jet 
terminates at the position of strong H2 bow-shocks, and 
shows a Hubble law velocity relation. Low-velocity CO 
traces a cavity that is very precisely located in the wake of 
the shocks. These observations strongly suggest that the 
propagation of one or several shocks in a protostellar jet 
entrain the ambient molecular gas and produces the low- 
velocity molecular outflow (see Sec. 3). With an estimated 
dynamical timescale of ~10^ years, HH211 is obviously 
an extremely young object. Other examples of such highly- 
colHmated, high-velocity jets include IRAS 041 66H-2706 
{Tafalla et al, 2004) and HH 212 {Lee et al, 2000) - fliese 
sources are or will be in the near future the subject of more 
detailed investigations. 

In at least SVS 13B {Bachiller et al, 1998, 2000), and 
NGC1333 IRAS 2 {J0rgensen et al, 2004), NGC1333 
IRAS 4 {Choi 2005) and HH211 {Chandler and Richer, 
2001; Hirano et al, 2006; Palau et al, 2006) the SiO 
emission traces the molecular jet and not the strong termi- 
nal shocks against the interstellar medium. This came as a 
surprise, as it seems to contradict the widely accepted idea 
that SiO is a tracer of outflow shocks, where the density 
is increased by several order of magnitudes (e.g., Martin- 
Pintado et al, 1992; Schilke et al, 1997; Gibb et al, 2004). 
The lack of significant SiO emission in the terminal shocks 
suggests that the formation process of this molecule has a 
strong dependence on the shock conditions (velocity, den- 
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sity) and/or outflow age (see Sec. 4.2). 

The exact nature of these CO and SiO molecular jets 
is not yet clear. Three basic scenarios could be invoked, 
in which the high- velocity CO and SiO molecules faj be- 
long to the actual protostellar jet, (b) are entrained along 
the jet in a turbulent cocoon (e.g., Stahler, 1994; Raga et 
ah, 1995), or (c) are formed/excited in shocks that are prop- 
agating down the jet ("internal working surfaces", Raga and 
Cabrit, 1993). This latter scenario would reconcile the ob- 
servation of SiO in the jet and the shock-tracer nature of 
this molecule. The predictions of these three cases, both 
in terms of line properties and observed morphologies, are 
somewhat different but the current observations have not 
yet provided a clear preference for one of these scenarios. 

2.7.2. More complex structures. Not all sources have 
structures as simple or unperturbed as the molecular jets 
discussed above. CO observations have also revealed a 
number of more complex outflow properties. 

Episodic ejection events seem to be a common prop- 
erty of young molecular outflows. In sources such as, 
e.g., LI 157 {Gueth et al, 1998) and IRAS 04239+2436 
(HH300, Arce and Goodman, 2001b), a limited number 
(2 to 5) of strong ejection events have taken place, each of 
them resulting in the propagation of a large shock. Morpho- 
logically, the flow is therefore the superposition of several 
shocked/outflowing gas structures, while position-velocity 
diagrams show multiple "Hubble wedges" (i.e., a jagged 
profile; Arce and Goodman, 2001a). In most of the sources, 
if several strong shocks are not present, a main ejection 
event followed by several smaller, weaker shocked areas 
are observed (e.g., L1448: Bachiller et al, 1990; HH 111: 
Cernicharo et al, 1999; several sources: Lee et al, 2000, 
2002). As noted before, even the molecular jets could in- 
clude several internal shocks. Altogether, these properties 
suggest that the ejection phenomenon in young outflows is 
intrinsically episodic, or — a somewhat more attractive pos- 
sibility — could be continuous but include frequent ejection 
bursts. This could be explained by sudden variations in the 
accretion rate onto the forming star, that result in variations 
of the velocity of the ejected matter, hence the creation of a 
series of shocks. 

Precession of the ejection direction has been estabhshed 
in a few sources, like Cep E {Eisloffel et al, 1996), and 
L 1 157 {Gueth et al, 1996, 1998). In several other objects, 
the observations reveal bending or misahgnment between 
the structures within the outflows (see e.g., Lee etal, 2000, 
2002). In fact, when observed at the angular resolution 
provided by millimeter interferometers, many well-defined, 
regular bipolar outflows mapped with single-dish telescopes 
often reveal much more complex and irregular structures, 
which indicate both temporal and spatial variations of the 
ejection phenomenon. 

Quadrupolar outflows are sources in which four lobes 
are observed, and seem to be driven by the same proto- 
stellar condensation. Several scenarios were proposed to 
explain these peculiar objects: two independent outflows 
(e.g., Anglada et al, 1991; Walker et a/., 1993); one sin- 



gle flow with strong limb-brightening, which would thus 
mimic four lobes (e.g., Avery et al, 1990); a single outflow 
but with a strong precession of the ejection direction (e.g., 
Ladd and Hodapp, 1997). The angular resolution provided 
by recent interferometric observations have clearly favored 
the first hypothesis in at least two objects (HH288, Gueth 
et al, 2001 ; L 723, Lee et al, 2002). In botii cases, the two 
outflows are driven by two independent, nearby protostars, 
located in the same molecular core. It is however unclear 
whether the sources are gravitationally bound or not. 

2.7.5. Time evolution. There is increasing evidence 
that outflow coUimation and morphology changes with time 
(e.g., Lee et al, 2002; Arce and Sargent, in preparation). 
The youngest outflows are highly collimated or include 
a very colhmated component, strongly suggesting that jet 
bow shock-driven models are appropriate to explain these 
objects. Older sources present much lower coUimation fac- 
tors, or — a somewhat more relevant parameter — wider 
opening angles, pointing towards wide-angle, wind-driven 
outflows (see Sec. 3.1.1). In fact, neither the jet-driven nor 
the wind-driven models can explain the range of morpho- 
logical and kinematic properties that are observed in all out- 
flows (see Sec. 3.2). This was noted by Cabrit et al. (1997), 
who compared outflow observations to morphologies and 
PV diagrams predicted by various hydrodynamical models. 
More recently, a similar conclusion was obtained by Lee et 
al (2000, 2001, 2002) from interferometric observations 
of 10 outflows. One attractive scenario to reconcile all ob- 
servations is to invoke the superposition of both a jet and a 
wind component in the underlying protostellar wind and a 
variation in time of the relative weight between these two 
components. One possible explanation for this scenario is 
that at very early ages only the dense collimated part of the 
wind can break out of the surrounding dense infalling en- 
velope. As the envelope loses mass, through infall and out- 
flow entrainment along the axis (see Sec 4. 1), the less dense 
and wider wind component will break through, entraining 
the gas unaffected by the collimated component, and will 
eventually become the main component responsible for the 
observed molecular outflow. 

2.2 Outflows from high-mass protostars 

Outflows from more luminous protostars have received 
increasing attention in recent years with the result that we 
now have a more consistent understanding of massive out- 
flow properties and their relationship to outflows from lower 
luminosity objects (see, e.g., recent reviews by Churchwell, 
1999; Shepherd, 2003, 2005; and Cesaroni, 2005). 

Outflows from mid- to early-B type stars have mass out- 
flow rates 10~^ to a few x 10"'^ M© yr~^, momentum rates 
10~* to 10~^ M0kms~^yr~^, and mechanical luminos- 
ity of 10~^ to 10^ L0. O stars with bolometric luminosity 
(Lfeoi) of more than Lq generate powerful winds with 
wind opening angle of about 90° within 50 AU of the star 
(measured from water masers in and along the flow bound- 
aries and models derived from ionized gas emission ob- 
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served with resolutions of 20-100 AU, e.g., Orion: Green- 
hill et al., 1998; MWC349A: Tafoya et al., 2004). The 
accompanying molecular flows can have an opening angle 
of more than 90° (measured from CO outflow boundaries 
1000 AU to 0.1 pc from the protostar). The flow momen- 
tum rate (> 10~^ kms~^yr~^) is more than an order 
of magnitude higher than what can be produced by steUar 
winds and the mechanical luminosity exceeds 10^ Lq (e.g., 
Churchwell, 1999; Gamy and Lizano, 1999). 

Outflows from early-B and late O stars can be well- 
collimated (collimation factors greater than 5) when the dy- 
namical times scale is less than ~ 10"^ years. For a few early 
B (proto)stars with outflows that have a well-defined jet, the 
jet appears to have adequate momentum to power the larger 
scale CO flow, although this relation is not as well estab- 
lished as it is for lower luminosity sources. For example, 
IRAS 20126-1-4104 has a momentum rate in the SiO jet of 

2>< 10"^ ( sfo/H" ) MQkms^Vi""^ while the CO momen- 
tum rate is 6 x 1O~^M0 km s~^yr~^ (Cesawni et al, 1999; 
Shepherd et al., 2000). Although the calculated momentum 
rate in the SiO jet is adequate to power the CO flow, the 
uncertainties in the assumed SiO abundance makes this dif- 
ficult to prove. Another example is IRAS 18151-1208 in 
which the H2 jet appears to have adequate momentum to 
power the observed CO flow (Beuther et al., 2002a; Davis 
et al., 2004). A counter example may be the Ceph A HW2 
outflow because the momentum rate in the HCO+ outflow is 
20 times larger than that of the observed ionized jet. How- 
ever, the jet could be largely neutral or there may be an un- 
detected wide-angle wind component (Gomez et al., 1999). 

Wu et al. (2004) find that the average collimation factor 
for outflows from sources with L^oi > 10^ L© is 2.05 com- 
pared with 2.81 for flows from lower luminosity sources. 
This is true even for sources in which the angular size of 
the flow is at least five times the resolution. Table 1 of 
Beuther and Shepherd (2005) summarizes our current un- 
derstanding of massive outflows from low-spatial resolution 
single-dish studies and gives a summary of and references 
to 15 massive flows that have been observed at higher spa- 
tial resolution using an interferometer. Here, we discuss a 
few of these sources that illustrate specific characteristics of 
massive outflows. 

2.2.7. Collimated flows. The youngest early-B proto- 
stars (~ 10^ years or less) can be jet-dominated and can 
have either well-collimated or poorly collimated molecular 
flows. In a few sources, jets tend to have opening angles, 
a, between 25° and 30° but they do not re-colUmate (e.g., 
IRAS 20126+4104: Cesaroni et al, 1999; MoscadeUi et 
al., 2005; or IRAS 16547-4247: Rodriguez et al., 2005a). 
Other sources appear to generate well-colhmated jets 
(a few degrees) that look like scaled up jets from low- 
luminosity protostars (e.g., IRAS 05358+3543: Beuther et 
al, 2002b). All these sources are ^10^ years old — they 
have not yet reached the main sequence. In at least one case 
jet activity has continued as long as 10^ years, although the 
associated molecular flow has a large opening angle and 



complex morphology (HH80-81: Yamashita et al, 1989; 
Marti etal, 1993). 

One possible collimated outflow event may have been 
traced to a young 05 (proto)star in flie G5. 89-0.39 UC HII 
region. The05 star has a small excess at 3.5//m and is along 
the axis of two H2 knots that appear to trace a N-S molecu- 
lar flow along the direction of the UC HII region expansion 
(Puga et al, 2005). The N-S molecular flow is unresolved 
so it is not clear that it is collimated even if the H2 knots 
appear to trace a collimated outflow event. Although still 
circumstantial, the evidence is mounting that the 05 star in 
G5.89 produced the N-S outflow and thus is forming via 
accretion {Shepherd, 2005, and references therein). 

2.2.2. Poorly collimated flows. Poorly collimated molec- 
ular flows can be due to: 1) extreme precession of the jet as 
in IRAS 20126+4124 {Shepherd et al, 2000); 2) a wide- 
angle wind associated with a jet as in HH 80-81 {Yamashita 
et al, 1989) or perhaps Ceph A HW2 (e.g., Gomez et al, 
1999; Rodriguez etal, 2001); 3) a strong wide-angle wind 
that has no accompanying jet; or 4) an explosive event as 
seen in Orion {McCaughrean and Mac Low, 1997). In 
massive flows, collimation factors as high as 4 or 5 in the 
molecular gas can still be consistent with being produced 
by wind-blown bubbles if the cloud core is very dense and 
it is easier for the flow to break out of the cloud rather than 
widen the flow cavity. Once the flow has escaped the cloud 
core, the bulk of the momentum is ttansfered to the inter- 
clump medium. 

In at least some young early-B stars, both the ionized 
wind near the central source and the larger scale molec- 
ular flow are poorly collimated and there is no evidence 
for a well-collimated jet. Examples of sources that do not 
appear to have a colhmated jet powering the flow include 
G192.16-3.82 {Shepherd and Kurtz, 1999, and references 
therein), W75N VLA2 {Torrelles et al, 2003, and refer- 
ences therein), AFGL490 {Schreyer et al, 2006, and refer- 
ences therein) and the SiO flow in G5. 89-0. 39 (not related 
to the 05 star discussed above; Sollins et al, 2004; Puga 
et al, 2005). Sources with poorly collimated flows, no ev- 
idence for a jet and a good determination of the dynamical 
age show that the ages tend to be a few x 10^ years old and 
a UC HII region exists around a new ZAMS star. 

To date, extremely colhmated molecular outflows have 
not been observed toward sources earlier than BO. It is pos- 
sible that this is simply a selection effect because O stars 
form in dense clusters and reach the ZAMS in only a few 
X 10^ years. Thus, any collimated outflows may be con- 
fused by other flows. In a few cases, outflows appear to 
be due to a sudden explosive event such as that seen in 
Spitzer images of shocked gas in G34.26+0.15 {Church- 
well, personal communication) or the H2 fingers of Orion. 
There is now good evidence that Source I in Orion and the 
Becklin-Neugebauer (BN) object were within a few hun- 
dred AU from each other about 500 years ago {Rodriguez 
et al, 2005b). Such close encounters could disrupt the ac- 
cretion process and create an explosive outflow as seen in 
Orion (e.g., Bonnell et al, 2003). 
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2.2.3. Evolution. Early-B stars (Lboi ^ 10'* L0) gen- 
erate UC HII regions and reach the ZAMS in 5 - 9 x 10'' 
years while still accreting and generating strong molecu- 
lar outflows (e.g., Churchwell, 1999; Garay and Lizano, 
1999, and references therein). The duration of the accre- 
tion phase is about the same as in low-luminosity sources 
(e.g., 5 — 10 X 10^ years) yet the development of an HII 
region that expands to encompass the accretion disk mid- 
way through the formation process suggests that there is a 
sharp transition in the physical conditions at the base of the 
flow where material is lifted off the surface of the disk and 
collimated. 

Well-collimated molecular flows from massive proto- 
stars tend to be in systems with ages less than a few 
times 10"* years old where the central object has not yet 
reached the main sequence (e.g., IRAS 05358+3543 is 
well-collimated over approximately 1 pc). In these young 
sources the effects of increased irradiation on the disk and 
disk-wind due to the stellar radiation field are minimal. 
Poorly collimated flows (opening angle greater than 50° 
that show no evidence for a more collimated component) 
are associated with more evolved sources that have de- 
tectable UC HII regions and the central star has reached the 
main sequence. 

To account for the differences seen in flow morphologies 
from early B to late O stars Beuther and Shepherd (2005) 
proposed two possible evolutionary sequences which could 
result in similar observable outflow signatures. In Fig. 1 we 
show a schematic of the proposed sequences and explain 
how the observed outflow morphologies can be related to O 
and B star evolution. 



HMPO 



HC HII 



UC HII 




Fig . 1 . — Sketch of the proposed evolutionary outflow scenario put forth 
by Beuther and Shepherd (2005). The three outflow morphologies can be 
caused by two evolutionary sequences: (top) the evolution of a typical 
B I -type star from a high-mass protostellar object (HMPO) via a hyper- 
compact Hll (HC HII) region to an ultra-compact Hll (UC HII) region, and 
(bottom) the evolution of an 05-type star which goes through Bl- and OS- 
type stages (only approximate labels) before reaching its final mass and 
stellar luminosity. This evolutionary sequence appears to qualitatively fit 
the observations, yet it must be tested against both theory and observations. 



Once a massive OB star reaches the main sequence, the 
increased radiation from the central star generates signifi- 
cant Lyman continuum photons and will likely ionize the 
outflowing gas even at large radii. Inherently lower colli- 
mation of the ionized wind due to increased radiation pres- 
sure is suggested by the hydrodynamic simulations of Yorke 
and Sonnhalter (2002). However the radiation pressure is 
still too low by a factor of 10 to 100 to produce significant 
changes in the collimation of the observed molecular flows 
{Richer ef a/., 2000). 

The larger photon flux will also increase the ionization 
degree in the molecular gas and produce shorter ion-neutral 
collisional timescales. Thus, in principle, this could im- 
prove the matter-field coupling, even aiding MHD collima- 
tion. However, other effects are likely to counteract this. In 
particular, if the plasma pressure exceeds the magnetic field 
pressure and ions are well-coupled to the field, then the out- 
flowing, ionized gas may be able to drag the magnetic field 
lines into a less collimated configuration (see, e.g., Kdnigl, 
1999; Shepherd etai, 2003). 

Turbulence could also contribute to the decollimation of 
molecular outflows from massive OB protostars. Increased 
turbulence in the disk and outflow is expected to weaken 
the conditions for ideal MHD and hence weaken the colli- 
mation effect. Turbulence could be due to higher accretion 
disk to stellar mass ratios {Mdisk > 0.37\f*) making disks 
susceptible to local gravitational instabilities, increased ra- 
diation pressure and high plasma temperatures. If the ions 
and neutrals are not well coupled in a turbulent flow then 
ideal MHD begins to break down and magnetic diffusiv- 
ity could significantly decollimate the molecular outflow 
(see, e.g., Fendt and Cemeljic, 2002). Further, simulations 
by Fendt and Cemeljic find that the toroidal magnetic field 
component, S^, decreases with increased turbulence. Since 
is the collimating magnetic component (e.g., Pudritz & 
Banerjee, 2005), such a decrease in may contribute to 
the lower observed collimation for more evolved massive 
molecular outflows. 

3. MOLECULAR OUTFLOW MODELS 
3.1 General Overview of Models 

Several outflow models have been proposed to explain 
how molecular outflows from protostars are formed. Cur- 
rently, outflow models can be separated into four broad 
classes (Cabrit et ai, 1997): (1) wind-driven shells, (2) jet- 
driven bow shocks, (3) jet-driven turbulent flows, and (4) 
circulation flows. In the first three, molecular outflows rep- 
resent ambient material that has been entrained by a wide- 
angle wind or accelerated by a highly collimated jet. In the 
last class of models, molecular outflows are produced by 
deflected infalling gas. Most of the work has concentrated 
on simulating outflows specifically from low-mass proto- 
stars, and little work has been done on modeling outflows 
from high-mass stars. Many flow properties, in particu- 
lar the CO spatial and velocity structure, are broadly simi- 
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lar across the entire luminosity range (Richer et ai, 2000), 
suggesting that similar mechanisms may be responsible for 
the production of molecular outflows from both low- and 
high-mass systems. Recent results from simulation work on 
the disk/outflow connection (Pudritz and Banarjee, 2005) 
as well as from observations {Zhang et al., 2002; Beuther 
et al, 2004) further indicate that molecular outflows from 
massive stars may be approximately modeled as scaled-up 
versions of their lower mass brethren. 

In the past, most studies used analytical models to try to 
explain the outflow morphology and kinematics. However, 
in the last decade computational power has increased suffi- 
ciently to allow for multidimensional hydrodynamical (HD) 
simulations of protostellar outflows that include a simple 
molecular chemical network. Numerical modeling of the 
molecular cooling and chemistry, as well as the hydrody- 
namics, is required in these systems, which are described by 
a set of hyperbohc differential equations with solutions that 
are usually mathematically chaotic and cannot be treated 
analytically. Treatment of the molecular cooling and chem- 
istry facihtates a comparison of the underlying flow with 
observational quantities (for example, the velocity distribu- 
tion of mass vs. CO intensity, the temperature distribution 
of the outflowing gas, and the H2 1-0 S(l) maps). 

3.1.1. Wind-driven shell models. In the wind-driven 
shell model, a wide-angle radial wind blows into the strati- 
fied surrounding ambient material, forming a thin swept-up 
shell that can be identified as the outflow shell {Shu et al, 
1991; Li and Shu, 1996; Matzner and McKee, 1999). In 
these models, the ambient material is often assumed to be 
toroidal with density pa = Pao sin^ while the wind is 
intrinsically stratified with density p,^ = (0^o/(r^ sin^ 6*), 
where pao is the ambient density at the equator and pwo is 
the wind density at the pole {Lee et al., 2001). This class of 
models is attractive as it particularly explains old outflows 
of large lateral extents and low colhmation. 

In recent years, there have been a few efforts to model 
wide angle winds numerically. Lee et al. (2001) performed 
numerical HD simulations of an atomic axisymmetric wind 
and compared it to simulations of bow shock-driven out- 
flows. Their wide-wind models yielded smaller values of 7 
(see Sec. 2) over a narrower range (1.3-1.8), as compared 
to the jet models (1.5-3.5). Raga et al. (2004b) have in- 
cluded both wide angle winds and bow shock models in a 
study aimed at reproducing features of the southwest lobe 
of HH 46/47, with the result that a jet model is able to match 
enough features that they feel that it is not necessary to in- 
voke a wide angle wind (although it produces a reasonable 
fit to the observations). In simulations by Delamarter et al. 
(2000) the wind is assumed to be spherical, even though 
the physical origin of such a wind is not yet clear, and it is 
focused towards the polar axis by the density gradients in 
the surrounding (infalling) torus-like environment. In these 
models the low- velocity 7 ranges from approximately 1.3 
to 1.5, similar to other studies. The MHD simulations per- 
formed by Gardiner et al. (2003) show that winds that have 
a wide opening angle at the base can produce a dense jet- 



Uke structure downstream due to MHD coUimation. Very 
recently, axisymmetric winds have been modeled with a 
code that includes molecular chemistry and cooling as well 
as Adaptive Mesh Refinement (AMR) {Cunningham et al, 
2005). These last two studies produce satisfactory gen- 
eral outflow lobe appearance, however, no mass-velocity, 
position-velocity maps, or channel maps have been gener- 
ated to compare with observations. 

3.1.2. Turbulent jet model. In the jet-driven tur- 
bulent model, Kelvin-Helmholtz instabilities along the 
jet/environmental boundary lead to the formation of a turbu- 
lent viscous mixing layer, through which the cloud molec- 
ular gas is entrained {Canto and Raga, 1991; Raga et al, 
1993; Stahler, 1994; Lizano and Giovanardi, 1995; Canto 
et al, 2003, and references therein). The mixing layer 
grows both into the environment and into the jet, and even- 
tually the whole flow becomes turbulent. Discussion of the 
few existing numerical studies that investigate how molec- 
ular outflows are created by a turbulent jet is presented 
in a recent review by Raga et al (2004a), who cite the 
"Torino group" as the only simulations with predictions for 
atomic (e.g.. Ha, [SII]) emission {Micono et al, 1998). 
The radiatively cooled jet simulations reproduce the bro- 
ken power law behavior of the observationally determined 
mass-velocity distribution, even though molecular chem- 
istry or cooling is not included {Micono et al, 2000). How- 
ever, these models produce decreasing molecular outflow 
momentum and velocity with distance from the powering 
source — opposite to that observed in most molecular out- 
flows. An analytical model using Kelvin-Helmhotz insta- 
bilities has recently been proposed by Watson et al (2004) 
to explain entrainment of cloud material by outflows from 
high-mass stars. 

3.1.3. Jet bow shock model. In the jet-driven bow shock 
model, a highly collimated jet propagates into the surround- 
ing ambient material, producing a thin outflow shell around 
the jet {Raga and Cabrit, 1993; Masson and Chemin, 
1993). The physical origin of the jet is currently unclear 
and could even be considered as an extreme case of a highly 
colUmated wide-angle wind without a tenuous wide-angle 
component. As the jet impacts the ambient material, a pair 
of shocks, a jet shock and a bow shock, are formed at the 
head of the jet. High pressure gas between the shocks is 
ejected sideways out of the jet beam, which then interacts 
with unperturbed ambient gas through a broader bow shock 
surface, producing an outflow shell surrounding the jet. An 
episodic variation in the mass-loss rate produces a chain of 
knotty shocks and bow shocks along the jet axis within the 
outflow shell. Recent analytical models without magnetic 
field include Wilkin (1996), Zhang and Zheng (1997), Smith 
etal. {\991), O striker et al. {2001), and Downes and Cabrit 
(2003). 

There have been two recent sets of efforts (by two differ- 
ent groups) to model molecular protostellar jets numerically 
in two or three spatial dimensions, where the mass- velocity 
and position-velocity have routinely been measured. In 
these simulations, a tracer associated with molecular hy- 
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drogen is followed. However, each group approaches this 
problem in a different way, with each approach having 
its own advantages and disadvantages. In an effort to re- 
solve the post-shock region, Downes and Ray (1999), and 
Dowries and Cabrit (2003) have simulated relatively low 
density, axisymmetric (two-dimensional) fast jets. Alter- 
natively, recognizing that observed flows associated with 
Class sources have a higher density and a complex appear- 
ance. Smith and Rosen have extended the work of Suttner 
etal. (1997) and Volker et al. (1999) by further investigat- 
ing sets of fully three-dimensional flows (e.g., Rosen and 
Smith, 2004a). The main disadvantage of this approach is 
that with such high densities the post-shock region will nec- 
essarily be under-resolved, especially in three-dimensional 
flows. Both the Downes and Smith groups have included 
molecular hydrogen dissociation and reformation as well 
as ro-vibrational cooling in their hydrodynamical simula- 
tions, although the treatment of this coohng is quite differ- 
ent in each group. One example is that the Downes group 
turns off all cooling and chemistry below 1000 K, while the 
Smith and Rosen simulations (explained in detail in Smith 
and Rosen, 2003) include coohng and chemistry calcula- 
tions at essentially all temperatures (albeit with an equilib- 
rium assumption for some reactions). The jet flows them- 
selves enter the grid from a limited number of zones at one 
side of the computational domain, with densities and tem- 
peratures that are constant radially (a top hat profile) and 
over time. Both groups usually model the jet as nearly com- 
pletely molecular — even though there are arguments sug- 
gesting that the jet will not initially be molecular, and that 
H2 might subsequently form on the internal working sur- 
faces of the jet (Raga et al, 2005). The initial jet velocities 
of the Downes, and Smith and Rosen groups are varied with 
shear, pulsation, and, in the three dimensional simulations, 
with precession. 

These different approaches have yielded different slopes 
for the computed CO intensity-velocity plots. The Downes 
group results have tended to be steeper and closer to the 
nominal value of 7 = 2, while the standard Rosen and Smith 
case has a value near 1. Much of this difference can be at- 
tributed to the difference in jet-to-ambient density ratio (see 
Rosen and Smith, 2004a), which is 1 in the Downes standard 
case, and 10 in the Rosen and Smith standard case. The 
value of 7 has been shown in these simulations to evolve 
over time, with steeper slopes associated with older flows. 
Most of these simulations are quite young, but there has 
been a recent effort to run the simulations out to t = 2300 
yr (Keegan and Downes, 2005). They confirm the steepen- 
ing of the mass- velocity slope up to t = 1600 yr (when 7 
= 1.6), and then it becomes roughly constant. The Smith 
and Rosen group have investigated whether fast (Rosen and 
Smith, 2004b) or slow (Smith and Rosen, 2005) precession 
has an effect on the mass-velocity slopes. While the sim- 
ulations with fast precessing jets show a dependence of 7 
on the precession angle (generally increasing 7 with the an- 
gle), some of this dependence was reduced in the slowly 
precessing cases. However, at this time only very young (t 



< 500 yr) precessing sources have been simulated. 

The initially molecular jet simulations that include pe- 
riodic velocity pulses exhibit position-velocity plots with 
a sequence of Hubble wedges, similar to that observed in 
molecular outflows produced by an episodic protostellar 
wind (see Sec. 2.1.2). Where computed, velocity channel 
maps in CO from molecular jet simulations, as in Rosen 
and Smith (2004a), have a morphology similar to that of 
many sources (e.g., HH 211, Gueth and Guilloteau, 1999), 
i.e. reveahng the knots within the jet at high velocities and 
showing the overaU shape of the bow shock at low veloci- 
ties. 

Some recent studies show the need to expand the in- 
terpretation of molecular outflow observations beyond the 
simulated H2 and CO emission from the numerical mod- 
els discussed above. For example, the work of Lesaffre et 
al. (2004) includes more complex chemistry in one dimen- 
sion, focusing on the unstable nature of combined C and J 
shocks. Also, radiation transfer with a complex chemistry 
has been simulated for a steady three dimensional (jet) flow, 
with a focus on HCO+ emission (Rawlings et al., 2004). 

In addition, magnetic field effects have been included 
in atomic protostellar jets that are axisymmetric (Gardiner 
et at., 2000; Stone and Hardee, 2000) and fully three- 
dimensional (Cerqueira and de Gouveia dal Pino, 1999, 
2001) and even molecular axisymmetric protostellar jets 
(O' Sullivan and Ray, 2000). These studies show significant 
differences compared to simulations of jets without mag- 
netic fields. For example, magnetic tension, either along the 
jet axis or as a hoop stress from a toroidal field, can help col- 
limate and stabihze the jet — ^though some of the additional 
stability is mitigated in a pulsed jet. Some of the differences 
between pure HD and MHD simulations that show up in the 
axisymmetric cases are less prominent in three dimensional 
simulations (Cerqueira andde Gouveia dal Pino, 2001). 

3.1.4. Circulation models. In circulation models the 
molecular outflow is not entrained by an underlying wind 
or jet, it is rather formed by infalling matter that is deflected 
away from the protostar in a central torus of high MHD 
pressure through a quadrupolar circulation pattern around 
the protostar, and accelerated above escape speeds by lo- 
cal heating (Fiege and Henriksen, 1996a,b). The molecu- 
lar outflow may still be affected by entrainment from the 
wind or jet, but this would be limited to the polar regions 
and it would not be the dominant factor for its acceleration 
(Levy et al, 1999, 2002). Circulation models may provide 
a means of injecting added mass into outflows from O stars 
where it appears unlikely that direct entrainment can supply 
all the observed mass in the flow (Churchwell, 1999). 

The most recent numerical studies of the circulation 
model have focused on a steady-state axisymmetric case, 
usually involving radiative heating, magnetic fields and 
Poynting flux (Lery, 2003). The addition of the Poynting 
flux in recent versions of this model has alleviated one of 
its major flaws (Lery et al., 2002), i.e. the inability in ear- 
lier models to generate an outflow of sufficient speed. The 
toroidal magnetic field in what is currently being called the 
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"steady-state transit model" assists in the formation of a col- 
limated fast moving flow {Combet et ai, 2006). 

3.2 Comparing Observations and Models 

In the past ten years, molecular outflows have been 
mapped at high angular resolutions with millimeter interfer- 
ometers, allowing us to confront the outflow models in great 
detail. A schematic of the predicted properties of molecular 
outflows produced by the different models discussed above 
is presented in Fig. 2. High-resolution molecular outflow 
observations can be used to compare the data with the out- 
flow characteristics shown in Fig. 2 in order to establish 
what model best fits the observed outflow. 

Molecular outflow properties predicted by different models 

Predicted property of molecular outflow along axis 



Model Wind Morphology Velocity Temperature Momentum" 




Assuming an underlying density distribution of r to r 



Fig. 2. — Observable molecular outflow properties predicted by the four 
leading broad classes of models: 1) turbulent jet (Canto and Raga, 1991; 
Chernin and Masson, 1995; Bence el al, 1996); 2) jet bow shock (Chemin 
and Masson, 1995; CUjfe et a/., 1996; Hatchell et al, 1999; Lee et al, 
2001); 3) wide-angle wind (Li and Shu, 1996; Lee et al, 2001); and 4) 
circulation models (Fiege and Henriksen, 1996b; Levy et al, 1999). In the 
jet-driven bow shock model, an episodic variation in jet velocity produces 
an internal bow shock driving an internal shell, in addition to the termi- 
nal shock. This episodic variation can also be present in the other wind 
models, but in this figure the effects of an episodic wind are only shown 
for the jet bow shock model. This figure is based on Figure 1 of Arce and 
Goodman (2002b). 

Here we focus our attention on comparing observations 
with the jet-driven bow shock and wide-angle wind-driven 
models, as most of the numerical simulations concentrate 
on these two models and they are the most promising mod- 
els thus far. The predicted mass-velocity relationships in jet 
bow shock and wide-angle wind models have a slope (7) 
of 1-4, in tune with observations. Each model predicts a 
somewhat different position-velocity (PV) relation that can 
be used to differentiate between these two leading molecu- 
lar outflow driving mechanisms (Cabrit et al, 1997; Lee et 
al, 2000, 2001). 

3.2.1. Jet-driven bow shock models vs. observations. 
Current jet-driven bow shock models can qualitatively ac- 
count for the PV spur structure (where the outflow veloc- 
ity increases rapidly toward the position of the internal and 



leading bow shocks, see Fig. 2), the broad range of CO 
velocities near H2 shocks, and the morphological relation 
between the CO and H2 emission seen in young and colli- 
mated outflows. These models are able to produce the ob- 
served outflow width for highly coUimated outflows, such 
as L1448, HH211 and HH212 (Bachiller et al, 1995; 
Gueth and Guilloteau, 1999; Lee et al, 2001). However, 
jet-driven bow shock models have difficulty producing the 
observed width of poorly collimated outflows, like RNO 9 1 , 
VLA 05487, and L 1221 (Lee et al, 2000, 2001, 2002). Jet 
models produce narrow molecular outflows mainly because 
the shocked gas in the bow shock working surfaces lim- 
its the transverse momentum (perpendicular to the jet-axis) 
that can be delivered to the ambient medium. In numerical 
simulations of jets, the width of the outflow shell is mainly 
determined by the effects of the leading bow shock from the 
jet's first impact into the ambient material (e.g., Suttner et 
al, 1997; Downes and Ray, 1999; Lee et al, 2001). While 
the jet penetration into the cloud increases roughly linearly 
with time, the width only grows as the one-third power of 
time (Masson and Chernin, 1993; Wilkin, 1996; Ostriker et 
al, 2001). 

Jets also have difficulty producing the observed outflow 
momenta. The transverse momentum of the outflow shell is 
acquired primarily near the jet head where the pressure gra- 
dient is large, and the mean transverse velocity of the shell, 
vh, can be approximated by vr ~ l3cs{R^j / R^), where R 
and Rj are the outflow and jet radius, respectively, and ficg 
is the velocity of the gas ejected from the working surface 
(Ostriker et al, 2001). For example, in a 10,000 AU-wide 
molecular outflow driven by a 150 AU jet, and assuming 
(ics = 32 km s~^, the expected mean transverse veloc- 
ity of the shell is only 0.03 km s~^. As a result, if out- 
flows were driven by a steady jet, the wide portions of out- 
flow shells would exhibit extremely low velocities and very 
small momenta. This is inconsistent with the observations, 
especially in the wider flows where the well-defined cavity 
walls have appreciable velocities (e.g., B5-1RS1: Velusamy 
andLanger, 1998; RN091: Lee et al, 2002; L1228: Arce 
and Sargent, 2004). 

Systematic wandering of the jet flow axis has been ar- 
gued to occur in several outflows based on outflow mor- 
phology, e.g., IRAS 20126H-4104 (Shepherd et al, 2000) 
and L 1157 (Bachiller et al, 2001). This may mitigate 
the above discrepancies. The width and momentum of the 
outflow shell can increase because a wandering jet has a 
larger "effective radius" of interaction and can impact the 
outflow shell more directly (Raga et al, 1993; Cliffe et al, 
1996). Some simulations show hints of widening by jet 
wandering (Vdlker et al, 1999; Rosen and Smith, 2004a; 
Smith and Rosen, 2005), but some show that a wandering 
jet could produce a smaller width than a steady jet (Raga et 
al, 2004b). Further calculations are needed to ascertain if 
motion of the jet axis at realistic levels can improve quanti- 
tative agreement with observed outflow features. 

3.2.2. Wide-angle wind models vs. observations. Wide- 
angle winds can readily produce CO outflows with large 
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widths but have trouble producing other commonly ob- 
served features. In this model, the outflow velocity also 
increases with the distance from the source, showing a lobe 
PV structure tilted with inclination that exhibits only a small 
velocity range at the tip. If the tip is not observed, the PV 
structure appears as a tilted parabola (see Fig. 2). As dis- 
cussed in Sec. 3.1.1, most wind-driven models, assume 
the protostellar wind density depends on the angle from 
the pole (9). If the wind velocity has a small, or no, de- 
pendence on 6, and assuming a density stratification sim- 
ilar to that proposed by Li and Shu (1996), then the out- 
flow width, W, can be expressed in terms of the ratio of 
wind to ambient density at the equator, [pwo/Pao), the 
wind velocity at the pole, v^o, and the outflow age, t, as 
W « {pwo/ PaoY^'^Vyjot {Lee etal., 2001). For (pwo/Pao) 
between 10~^ and 10~^, a 100 km s~^ wind can produce 
an outflow width of 0.1 to 0.2 pc in 10'* years. Thus, the 
wind-driven model can produce widths consistent with ob- 
served molecular outflows in about 10* years. However, 
these models have problems producing discrete bow shock 
type features in the entrained molecular gas, as seen in 
many high-resolution maps of CO outflows (e.g., Lee et al, 
2000, 2002), and discrete position-velocity spur structures 
(and Hubble wedges). These features are hard to generate 
as the wide wind impacts all locations on the shell. Models 
of wide-angle pulsed winds produce a series of flat inter- 
nal shocks within the outflow shell {Lee et al, 2001), in- 
consistent with the curved internal H2 bow shocks typically 
observed in episodic outflows (see Sec 2.1). 

One possible solution to these problems is to require the 
winds to have a coUimated core with a strong velocity gra- 
dient with respect to 9. A disk-wind driven from a large 
range of radii may have velocity strongly decreasing toward 
equatorial latitudes, because the asymptotic velocity on a 
given streamline in an MHD wind is characteristic of the 
Keplerian speed at the streamhne's footpoint (see chapter 
by Pudritz et al.). Further work is needed to study whether 
this sort of modification can produce the observed outflow 
features. 

3.2.3. A synthesis with an evolutionary scenario. A 
model which combines attributes of the jet and wide-angle 
wind models is arguably the best match to the available 
CO outflow data. A two component protostellar wind 
may be produced, for example, by a slow disk wind and 
a fast central disk-driven jet or X-wind (arising from the 
magnetosphere-disk boundary region) . The disk wind could 
help collimate the X-wind into the jet component (Ostriker, 
1997) and provide a slow wide-angle component that drives 
the outflow width and momentum (see chapter by Shang et 
at.). 

Observational support for the synthesis model exist 
at different wavelengths. There is mounting evidence 
from millimeter observations that the morphology of some 
molecular outflows is better explained with a "dual-wind" 
model (e.g., Yu et al, 1999; Arce and Goodman, 2002a; 
Arce and Sargent, 2004). In the optical, the forbidden 
emission line profiles of T Tauri stars show two velocity 



components: a high-velocity component that is argued to 
arise in a jet and a low- velocity component that might result 
from a disk- wind {Kwan and Tademaru, 1995; chapter by 
Ray et al). A possible scenario is that the main driving 
agent producing most of the observed molecular outflow 
may change over the time, as discussed in Sec. 2. Nu- 
merical simulations of an evolving dual-wind model will 
be critical to study whether this proposed scenario can re- 
produce the wide range of observed features in molecular 
outflows from low- and high-mass protostars. 

4. IMPACT OF OUTFLOWS ON SURROUNDING 
ENVIRONMENT 

4.1 Physical Impact 

Outflows from newborn stars inject momentum and en- 
ergy into the surrounding molecular cloud at distances rang- 
ing from a few AU to up to tens of parsecs away from the 
source. Historically, most studies have concentrated on the 
interaction between the outflow and the surrounding core 
(^ 0. 1 to 0.3 pc) as these scales can easily be observed with 
single-dish telescopes in the nearby (< 1 kpc ) star forming 
regions. More recently, studies using millimeter interfer- 
ometer array and single telescopes with focal-plane arrays 
have been crucial in the understanding of the outflow's im- 
pact at smaller (< 0. Ipc) and larger (> 1 pc) scales, respec- 
tively. 

4.7.7. Outflow-envelope interactions. Protostellar winds 
originate within a few AU of the star (see chapter by Ray 
et al), and so they are destined to interact with the dense 
circumstellar envelope — the primary mass reservoir of the 
forming star, with sizes in the range of 10^ to 10* AU. 
In fact, survey studies of the circumstellar gas within 10* 
AU of low-mass YSOs show outflows contribute signifi- 
cantly to the observed mass-loss of the surrounding dense 
gas (from about 1G~^ to 1G~* Moyr"*, depending on the 
protostar's age) and indicate there is an evolution in the 
outflow-envelope interaction (e.g.. Fuller and Ladd, 2002; 
Arce and Sargent, in preparation). As shown below, de- 
tailed studies of individual sources corroborate these re- 
sults. The powerful outflows from low-mass class sources 
are able to modify the distribution and kinematics of the 
dense gas surrounding a protostar, as evidenced in L 1157 
{Gueth et al, 1997; Beltrdn et al, 2004b), and RN043 
{Arce and Sargent, 2005) where molecular line maps show 
the circumstellar high-density gas has an elongated struc- 
ture and a velocity gradient, at scales of 4000 AU, along the 
outflow axis. Similarly, in IRAM0491 {Lee et al, 2005) 
and HH 212 {Wiseman et al, 2001) the dense gas traced by 
N2H+ and NH3, respectively, exhibit blue- and red-shifted 
protrusions extending along the blue and red outflow lobes, 
evidence that there are strong outflow-envelope interactions 
in these class sources. These results clearly show that, in- 
dependent of the original (i.e., pre-protostellar outflow) un- 
derlying circumstellar matter distribution, young outflows 
entrain dense envelope gas along the outflow axis. 
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Although not as powerful as those of class sources, the 
wide-angle outflows typically observed in class 1 sources 
(with opening angles of > 90°) are capable of constrain- 
ing the infalling envelope to a Umited volume outside the 
outflow lobes, as seen in the L 1228 (Arce and Sargent, 
2004) and B5-IRS1 (Velusamy and Langer, 1998) outflows. 
The LI 228 outflow is currently eroding the surrounding en- 
velope by accelerating high-density ambient gas along the 
outflow-envelope interface and has the potential to further 
widen the cavities, as the outflow ram pressure is about a 
factor of 4 higher than the infall ram pressure (Arce and 
Sargent, 2004). In RNO 91, a class II source, the outflow 
exhibits an even wider opening angle of 160° that is ex- 
panding, and decreasing the volume of the infall region (Lee 
and Ho, 2005). 

Widening of the outflow opening angle with age appears 
to be a general trend in low-mass protostars and there is 
ample evidence for erosion of the envelope due to outflow- 
envelope interactions (Velusamy and Langer, 1998; Arce 
and Sargent, 2004; Arce, 2004; Lee and Ho, 2005; Arce 
and Sargent, in preparation). Thus, it is clear that even if the 
pre-protostellar outflow circumstellar distribution of matter 
has a lower density along the polar regions (i.e., the out- 
flow axis) as suggested by different models (i.e., Hartmann 
et al., 1996; Li and Shu, 1996), outflow-envelope interac- 
tions will have an impact on the subsequent circumstellar 
density distribution, as they will help widen the cavity and 
constrain the infall region. It is tempting to extrapolate and 
suggest that as a young star evolves further its outflow will 
eventually become wide enough to end the infall process 
and disperse the circumstellar envelope altogether. 

4.7.2. Outflow-core interactions. Strong evidence ex- 
ists for the disruptive effects outflows have on their parent 
core — the dense gas within 0.1 to 0.3 pc of the young star. 
Direct evidence of outflow-core interaction comes from the 
detection of velocity shifts in the core's medium and high- 
density gas in the same sense, both in position and velocity, 
as the high-velocity (low-density) molecular outflow traced 
by ^^CO (e.g., Tafalla and Myers, 1997; Dobashi and Ue- 
hara, 2001; Takakuwa et al., 2003; Beltrdn et al., 2004a). 
The high opacity of the ^^CO lines hampers the ability to 
trace low-velocity molecular outflows in high-density re- 
gions. Therefore, other molecular species like ^^CO, CS, 
C^^^O, NH3, CH3OH, and C3H2 are used to trace the high- 
density gas pertorbed by the underlying protostellar wind. 
The average velocity shifts in the dense core gas are typ- 
ically lower than the average velocity of the molecular 
(^^CO ) outflow, consistent with a momentum-conserving 
outflow entrainment process. In addition to being able to 
produce systematic velocity shifts in the gas, outflows have 
been proposed to be a major source of the turbulence in the 
core (e.g., Myers etal, 1988; Fuller and Ladd, 2002; Zhang 
etal, 2005). 

Outflows can also reshape the structure of the star- 
forming core by sweeping and clearing the surrounding 
dense gas and producing density enhancements along the 
outflow axis. The clearing process is revealed by the pres- 



ence of nebular emission resulting from the scattering of 
photons, from the young star, off of cavity walls created by 
the outflow (e.g., Yamashita et al, 1989; Shepherd et al., 
1998; Yu etal, 1999), or depressions along the outflow axis 
in millimeter molecular hne maps of high density tracers 
(e.g., Moriarty-Schieven and Snell, 1988; White and Frid- 
lund, 1992; Tafalla etal, 1997). Outflow-induced density 
enhancements (and shock-heated dust) in the core may be 
revealed by the dust continuum emission (e.g., Gueth et 
al, 2003; Beuther et al, 2004; Sollins et al, 2004). A 
change in the outflow axis direction with time, as observed 
in many sources (see chapter by Bally et al.) will allow an 
outflow to interact with a substantial volume of the core and 
be more disruptive on the dense gas than outflows with a 
constant axis (e.g.. Shepherd et al, 2000; Arce and Good- 
man, 2002a). By accelerating and moving the surrounding 
dense gas, outflows can gravitationally unbind a signifi- 
cant amount of gas in the dense core thereby hmiting the 
star formation efficiency of the dense gas (see Matzner and 
McKee, 2000). 

The study of Fuente et al (2002) shows that outflows 
appear to be the dominant mechanism able to efficiently 
sweep out about 90% of the parent core by the end of the 
pre-main sequence phase of young intermediate-mass (Her- 
big Ae/Be) stars. In addition, outflows from low- and high- 
mass protostars have kinetic energies comparable to the 
gravitational binding energy of their parent core, suggesting 
outflows have the potential to disperse the entire core (e.g., 
Tafalla and Myers, 1997; Tafalla et al, 1997). We may 
even be observing the last stages of the outflow-core inter- 
action in G192.16, a massive (early B) young star, where 
the dense core gas is optically thin and clumpy, and the 
ammonia core is gravitationally unstable {Shepherd et al., 
2004). However, further systematic observations of a sta- 
tistical sample of outflow-harboring cores at different ages 
are needed in order to fully understand the details of the 
core dispersal mechanism and conclude whether outflows 
can disperse their entire parent core. 

Theoretical studies indicate that shocks from a proto- 
stellar wind impacting on a dense clump of gas (i.e., a 
pre-stellar core) along the outflow's path can trigger col- 
lapse and accelerate the infall process in the impacted core 
{Foster and Boss, 1996; Motoyama and Yoshida, 2003). 
Outflow-triggered star-formation has been suggested in 
only a handful of sources where the morphology and veloc- 
ity structure of the dense gas surrounding a young protostar 
appears to be affected by the outflow from a nearby YSO 
{Girart et al, 2001; Sandell and Knee, 2001; Yokogawa et 
al, 2003). 

4.L3. Outflow-cloud interactions far from the source. 
Giant outflows from young stars of all masses are common, 
and they can interact with the cloud gas at distances greater 
than 1 pc from their source (Reipurth et al, 1997; Stanke et 
al, 2000). Outflows from low-mass protostars are able to 
entrain 0.1 to 1 M solar masses of cloud material, acceler- 
ate and enhance the linewidth of the cloud gas {Hence et al, 
1996; Arce and Goodman, 2001b), and in some cases their 
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kinetic energy is comparable to (or larger than) the turbu- 
lent energy and gravitational binding energy of their par- 
ent cloud {Arce, 2003). The effects of giant outflows from 
intermediate- and high-mass YSOs on their surroundings 
can be much more damaging to their surrounding environ- 
ment. Studies of individual sources indicate that giant out- 
flows are able to entrain tens to hundreds of solar masses, 
induce parsec-scale velocity gradients in the cloud, produce 
dense massive sheUs of swept-up gas at large (> 0.5 pc) 
distances from the source, and even break the cloud apart 
{Fuente et al., 1998; Shepherd et al., 2000; Arce and Good- 
man, 2002a; Benedettini et al, 2004). The limited number 
of studies in this field suggest that a single giant outflow 
has the potential to have a disruptive effect on their parent 
molecular cloud (e.g., Arce, 2003). Clearly, additional ob- 
servations of giant outflows and their clouds are needed in 
order to quantify their disruptive potential. 

Most star formation appears in a clustered mode and so 
multiple outflows should be more disruptive on their cloud 
than a single star. Outflows from a group of young stars 
interact with a substantial volume of their parent cloud by 
sweeping up the gas and dust into shells (e.g., Davis et al, 
1999; Knee and Sandell, 2000), and may be a considerable, 
albeit not the major, source of energy for driving the su- 
personic turbulent motions inside clouds (Yu et al., 2000; 
Williams et al., 2003; Mac Low and Klessen, 2004). It has 
also been suggested that past outflow events from a group 
of stars may leave their imprint on the cloud in the form of 
numerous cavities (e.g.. Bally et al, 1999; Quillen et al, 
2005). Very limited (observational and theoretical) work on 
this topic exists, and further observations of star forming 
regions with different environments and at different evolu- 
tionary stages are essential to understand the role of out- 
flows in the gaseous environs of young steUar clusters. 

4.2 Shock chemistry 

The propagation of a supersonic protostellar wind 
through its surrounding medium happens primarily via 
shock waves. The rapid heating and compression of the 
region trigger different microscopic processes — such as 
molecular dissociation, endothermic reactions, ice sublima- 
tion, and dust grain disruption — which do not operate in 
the unperturbed gas. The time scales involved in the heating 
and in some of the "shock chemistry" processes are short (a 
few 10^ to 10* yr), so the shocked region rapidly acquires a 
chemical composition distinct from that of the quiescent un- 
perturbed medium. Given the short shock coohng times (~ 
10^ yr, Kaufman and Neufeld, 1996), some of these high- 
temperature chemical processes only operate at the initial 
stages, as the subsequent chemical evolution is dominated 
by low temperature processes. This chemical evolution, the 
gradual clearing of the outflow path, and the likely intrinsic 
weakening of the main accelerating agent, all together make 
the important signatures of the shock interaction (including 
some of the chemical anomalies) vanish as the protostellar 
object evolves. Chemical anomalies found in an outflow 



can therefore be considered as an indicator of the outflow 
age (e.g., Bachiller et al, 2001). 

The chemical impact of outflows are better studied in 
outflows around Class sources with favorable orienta- 
tion in the sky (i.e., high inclination with respect to the 
Une of sight). With less confusion than that found around 
massive outflows, the shocked regions of low-mass, high- 
collimation outflows (which often adopt the form of well- 
defined bows) are well separated spatially with respect to 
the quiescent gas. Detail studies of these "simple" regions 
can help disentangle the effects of outflow shocks from 
other shocks in more complex regions — like in circum- 
stellar disks, where one expects to find outflow shock ef- 
fects blended with those produced by shocks triggered by 
the collapsing envelope (e.g., Ceccarelli etai, 2000). 

Shocks in molecular gas can be of C-type or of J-type, 
depending on whether the hydrodynamical variables change 
continuously across the shock front (e.g., Draine and Mc- 
Kee, 1993). C-shocks are mediated by magnetic fields act- 
ing on ions that are weakly coupled with neutrals, they 
are slow, have maximum temperatures of about 2000-3000 
K, and are non-dissociative. J-shocks are typically faster, 
and can reach much higher temperatures. The critical ve- 
locity at which the change between C- and J-regime is 
produced depends on several parameters such as the pre- 
shock density (Le Bourlot et al, 2002) and the presence 
of charged grains {Flower and Pineau des Forets, 2003), 
and it typically ranges from ~ 20 up to ~ 50 km s~^. 
J-shocks may also occur at relatively low velocities when 
the transverse component of the magnetic field is small 
{Flower et al., 2003). Recent infrared observations of sev- 
eral lines of H2, CO, H2O, and OH, and of some crucial 
atomic fines, have made possible the estimate of tempera- 
ture and physical conditions in a relatively large sample of 
outflows. It follows that the interpretation of the data from 
most shocked regions require a combination of C- and J- 
shocks (see Noriega-Crespo, 2002; van Dishoeck, 2004, for 
comprehensive reviews). Such a combination of shocks can 
be obtained by the overlap of multiple outflow episodes as 
observed in several sources, and/or by the bow shock ge- 
ometry which could generate J-shocks at the apex of the 
bow together with C-shocks at the bow flanks {Nisini et al, 
2000; O'Connell et al, 2004, 2005). C-shocks are particu- 
larly efficient in triggering a distinct molecular chemistry in 
the region in which the molecules are preserved and heated 
to ~ 2000-3000 K. Moreover, molecules can also reform in 
J-shocked regions when the gas rapidly cools, or in warm 
layers around the hottest regions. The main processes ex- 
pected to dominate this shock chemistry were discussed by 
Richer et al (2000). 

Comprehensive chemical surveys have been carried out 
in two prototypical Class sources (L1157: Bachiller and 
Perez-Gutierrez, 1997; BHR71: Garay et al, 1998). More 
recent observations, including high-resolution molecular 
maps, have been made for a sample of sources, for ex- 
ample: LI 157 {Bachiller et al, 2001), NGC1333 IRAS2 
{J0rgensen et al, 2004), NGC1333 IRAS 4 {Choi et al. 
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2004), NGC2071 (Gamy et al, 2000), Cep-A (Codella et 
al., 2005). These observations have revealed that there are 
important differences in molecular abundances in different 
outflow regions. Such variations in the abundances may be 
linked to the time evolution of the chemistry (Bachiller et 
al., 2001) and may also be related to variations in the abun- 
dance of the atomic carbon {Jprgensen et al, 2004). 

SiO exhibits the most extreme enhancement factors 
(up to ~ 10®) with respect to the quiescent unperturbed 
medium. Such high enhancements are often found close to 
the heads (bowshocks), and along the axes, of some highly 
colUmated outflows (e.g., Dutrey et al., 1997, and refer- 
ences therein; Codella et al., 1999; Bachiller et al, 2001; 
Garay et al, 2002; J0rgensen et al, 2004; Patau et al, 
2006, and references therein). Sputtering of atomic Si from 
the dust grains is at the root of such high SiO abundances 
(Schilke et al, 1997), a process which requires shock ve- 
locities in excess of ~ 25 km s~^. Accordingly, the SiO 
lines usually present broad wings and, together with CO, 
the SiO emission usually reaches the highest terminal ve- 
locities among all molecular species. Moreover, recent 
observations of several outflows have revealed the presence 
of a narrow (< 1 km s~^) SiO line component (Lefloch 
et al., 1998; Codella et al, 1999; Jimenez-Serra et al, 
2004). The presence of SiO at low velocities is not well 
understood. Plausible explanations include that this is the 
signature of a shock precursor component {Jimenez-Serra 
et al, 2004, 2005) or that SiO is indeed produced at high 
velocities and subsequently slowed down in time scales of 
~ 10" yr {Codella et al, 1999). 

CH3OH and H2CO are also observed to be significantly 
overabundant in several outflows, enhanced by factors of 
about 100 {Bachiller et al, 2001; Garay et al, 2000; Garay 
et al, 2002; J0rgensen et al, 2004; Maret et al, 2005). 
These two species are likely evaporated directly from the 
icy dust mantles, and in many cases the terminal veloci- 
ties of their line profile wings are significantly lower than 
that of SiO, probably because CH3OH and H2CO do not 
survive at velocities as high as those required to form SiO 
{Garay etal, 2000). Thus, an enhancement of CH3OH and 
H2CO with no SiO may indicate the existence of a weak 
shock. On the other hand, after the passage of a strong 
shock, and once the abundances of CH3OH, H2CO and 
SiO are enhanced in the gas phase, one would expect the 
SiO molecules to re-incorporate to the grains while some 
molecules of CH3OH and H2CO remain in the gas-phase, 
as these two molecules are more volatile than SiO (their 
molecular depletion timescales are about a few 10^ yr for 
densities of ~ 10^ cm~^). In this scenario enhancement of 
CH3OH and H2CO most likely may mark a later stage in the 
shock evolution than that traced by high SiO abundances. 

In several outflows HCO+ high velocity emission is 
only prominent in regions of the outflow which are rela- 
tively close to the driving sources {Bachiller et al, 2001; 
J0rgensen et al, 2004). In such regions, the HCO"*" abun- 
dance can be enhanced by a factor of '^20. This behav- 
ior can be understood if the HCO"*" that was originally pro- 



duced through shock-induced chemistry (e.g., Rawlings et 
al, 2004) is destroyed by dissociative recombination or by 
reaction with the abundant molecules of H2O {Bergin et 
al, 1998). Once the abundance of the gaseous H2O de- 
creases due to freeze-out, the abundance of HCO"*" may 
increase. A rough anti-correlation between CH3OH and 
HCO+ {J0rgensen et al, 2004) seems to support these ar- 
guments. In other cases, HCO+ emission is observed at 
positions close to HH objects that can be relatively distant 
from the driving sources. In fact, together with NH3, HCO+ 
is expected to be enhanced in clumps within the molecular 
cloud by UV irradiation from bright HH objects {Viti and 
Williams, 1999), an effect observed near HH2 according to 
Girartetal (2002). Nevertheless, Girartetal. (2005) have 
recently found that UV irradiation alone is insufficient to 
explain the measured HCO+ enhancements and that strong 
heating (as that caused by a shock) is also needed. 

The chemistry of sulfur bearing species is of special in- 
terest as it has been proposed to be a potential tool to con- 
struct chemical clocks to date outflows (and hence their pro- 
tostellar driving sources). The scenario initially proposed 
by a number of models is that H2S is the main reservoir 
of S in grain mantles, although recent observations seem 
to indicate that OCS is more abundant on ices than H2S 
{Palumbo et al., 1997; van der Tak et al, 2003). Once 
H2S is ejected to the gas phase by the effect of shocks, its 
abundance will rapidly decrease after 10^ yr (e.g., Charn- 
ley, 1997) due to oxidation with O and OH, thereby pro- 
ducing SO (first) and SO2 (at a later time). Models and 
observations indicate that the SO/H2S and SO2/H2S ratios 
are particularly promising for obtaining the relative age of 
shocks in an outflow (C/zarnZey, 1991 ; Hatchell et al, 1998; 
Bachiller et al, 2001; Buckle and Fuller, 2003). On the 
other hand, recent models by Wakelam et al (2004) have 
shown that the chemistry of sulfur can be more complex 
than previously thought since — among other reasons — 
the abundances of the sulfur-bearing species critically de- 
pend on the gas excitation conditions, which in turn depend 
on the outflow velocity structure. Wakelam et al (2005) 
used the SO2/SO and the CS/SO ratios to constrain the age 
of the NGC1333 1RAS2 outflow to < 5 x lO^ yr. A re- 
cent study by Codella et al. (2005) confirms that the use 
of the SO/H2S and SO2/H2S ratios is subject to important 
uncertainties in many circumstances, and that other molec- 
ular ratios (e.g., CH3OH/H2CS, OCS/H2CS) can be used as 
more effective chemical clocks to date outflows. 

Recent work has revealed that chemical studies can be 
useful for the investigation of interstellar gas structure. For 
instance, Viti et al. (2004) have recently shown that, if the 
outflow chemistry is dominated by UV irradiation, clump- 
ing in the surrounding medium prior to the outflow passage 
is needed in order to reproduce the observed chemical abun- 
dances in some outflows. We stress, however, that this re- 
sult depends on the chemical modeling and that more work 
is needed before it can be generalized. 
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5. FUTURE WORK 

We discussed how the high angular resolution obser- 
vations have revealed general properties and evolutionary 
ttends in molecular outflows from low- and high-mass pro- 
tostars. However, these results rely on a limited number of 
outflows maps, thus making any statistical analysis some- 
what dangerous. A large sample of fully mapped outflows 
at different evolutionary stages, using millimeter interfer- 
ometers, is needed to soundly establish an empirical model 
of outflow evolution, and the outflow's physical and chem- 
ical impact on its surroundings. Also, detail mapping of 
many outflows will enable a thorough comparison with dif- 
ferent numerical outflow models in order to study the out- 
flow entrairmient process. 

Further progress in our understanding of outflows is ex- 
pected from current or planned instrument developments 
that aim at improving both the sensitivity and the angular 
resolution, while opening new frequency windows. The 
soon to be implemented improvements to the TRAM Plateau 
de Bure interferometer — which include longer baselines, 
wider frequency coverage, and better sensitivity — as well 
as the soon to be operational Combined Array for Research 
in Millimeter- wave Astronomy (C ARMA) will allow multi- 
Une large-scale mosaic maps with 1" resolution (or less), 
required to thoroughly study the outflow physical proper- 
ties (e.g., kinematics, temperature, densities), the entrain- 
ment process and the different chemical processes along 
the outflows' entire extent. In addition, large-scale mosaic 
maps of clouds with outflows will allow the study of the im- 
pact of many outflows on their parent cloud. The Atacama 
Large Milhmeter Array (ALMA), presumably operational 
by 2012, will have the ability to determine high-fidelity 
kinematics and morphologies of even the most distant out- 
flows in our Galaxy as well as flows in near-by galaxies. 
The superb (sub-arcsecond) angular resolution will be par- 
ticularly useful to study how outflows are ejected from ac- 
cretion disks, how molecular gas is entrained in the out- 
flow, and the interaction between the molecular jet/outflow 
and the environment very close to the protostar (i.e., the in- 
falhng envelope, and protoplanetary disk). The Expanded 
Very Large Array (EVLA), expected to be complete in 
2012, will be critical to image the wide-opening angle, ion- 
ized outflow close to the powering source, and will allow 
sensitive studies of re-ionization events in jets, H2O masers 
and SiO(l-O) in outflows. 

New submillimeter facilities and telescopes under con- 
struction will soon provide sensitive observations of high 
excitation lines, important for the study of outflow driving 
and entrainment, as well as shock-induced chemical pro- 
cesses. The recently dedicated Submillimeter Array (SMA) 
is the first instrument capable of studying the warm molecu- 
lar gas in the CO(6-5) line, at (sub)arcsecond resolution, al- 
lowing to trace the outflow components closer to the driving 
source and closer to the jet axis than previously possible. 
Furthermore, the large bandwidth of the correlator allows 
for simultaneous multi-Une observations crucial for study- 



ing the various shock chemistry processes in the outflow. 
Also, the Herschel Space Observatory (HSO) will measure 
the abundances of shock tracers of great interest, in partic- 
ular water, which cannot be observed from the ground. 

In the near future, greater computing power will make 
possible larger scale numerical simulations that take advan- 
tage of adaptive grids, better and more complex cooUng and 
chemistry functions, and the inclusion of radiative transfer 
and magnetic fields. Given the wealth of high-resolution 
data that will soon be available, numerical studies will need 
to compare the simulated outflows with observations in 
more detail, using the outflow density, kinematics, temper- 
ature and chemical structure. In addition, simulations that 
run for far longer times (^ IC* — 10"'' yr) than current mod- 
els (~ 10^ yr) are needed to study the outflow temporal 
behavior and evolution. Advances in computing, perhaps 
including GRID technology, may even allow a version of a 
virtual telescope, where both numerical modelers and ob- 
servers can find the best fit from a set of models for different 
sources. 
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